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National  Center  for  Atmospheric  Research,1  P.O.  Box  3000,  Boulder,  CO  80307 

Abstract.  Magnetized  electron  plasmas  have  been  used  to  study  the  dynamics  of 
two-dimensional  vortices  in  a cloud  of  background  vorticity.  Experiments  have  shown 
that  background  vorticity  can  calm  chaotic  vortex  motion,  and  cool  a system  of  intense 
vortices  into  a crystal  equilibrium.  Further  experiments  have  shown  that  weak  vortices 
tend  to  migrate  to  extrema  in  the  background  vorticity  distribution.  The  electron 
plasma  experiments  have  motivated  new  theories  on  two-dimensional  vortex  dynamics 
with  background  vorticity.  The  experiments  and  related  theories  are  summarized. 


Energetic,  roughly  two-dimensional  (2D)  vortices  abound  in  nature  [1],  Examples 
are  oceanic  eddies,  tropical  cyclones,  and  the  Jovian  hurricanes.  Remarkably,  mag- 
netized electron  plasmas  are  excellent  systems  for  studying  high  Reynolds-number 
2D  vortex  dynamics  experimentally  [2].  Recent  electron  plasma  experiments  have 
shed  light  on  how  2D  vortices  interact  with  the  environmental  flow,  i.e.,  with  back- 
ground vorticity.  This  article  reviews  the  experiments,  and  related  theories. 


ELECTRON  PLASMA  EXPERIMENTS 

Figure  1 is  a schematic  diagram  of  the  experimental  device  [2].  A long  column  of 
electrons  is  confined  in  a hollow  cylindrical  conductor.  The  electrons  are  trapped 
axially  by  negative  dc  voltages  on  two  end-cylinders,  and  radially  by  a uniform 
axial  magnetic  field  Bz.  The  apparatus  is  equipped  to  destructively  measure  the 
z-averaged  electron  density  n(r,0,t).  By  raising  the  voltage  of  one  end-cylinder 
(right),  the  electrons  are  released  onto  a phosphor  screen.  The  resulting  density 
image  is  recorded  with  a CCD  camera.  The  evolution  of  n is  observed  by  a sequence 
of  experiments,  in  which  the  the  initial  conditions  are  the  same,  but  the  release- 
times  vary. 

Under  proper  conditions  [2],  the  electron  density  n is  governed  approximately  by 
a set  of  two-dimensional  fluid  equations,  known  as  2D  drift-Poisson  flow: 
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FIGURE  1.  Confinement  device  (Penning-Malmberg  trap)  and  imaging  diagnostic,  used  to 
study  2D  vortex  dynamics  with  an  electron  plasma. 


UlX  C 

- + v-Vn  = 0,  v=—V—<j)Xz,  V2<f>  = Aiteri.  (1) 

at  B 

Here,  v(r,9,t)  is  the  E x B drift  velocity  field,  and  <f>(r,0,t)  is  a 2D  electrostatic 
potential.  The  boundary  condition  on  the  electrostatic  potential  is  (j)  = 0 at  the 
wall  of  the  confinement  cylinder. 

Defining  the  streamfunction,  xj)(r,0,t)  = cfy/B,  and  the  vortieity,  (( r,9,t ) = 
z • V x v , we  may  recast  2D  drift-Poisson  flow  as  follows: 

~+tT-V(  = 0,  v = -VV' x z,  V20  = £.  (2) 

By  amazing  coincidence,  Eqs.  (2)  are  the  2D  Euler  equations,  which  also  govern 
two-dimensional  flows  in  uniform-density  inviscid  incompressible  fluids. 

Comparing  (2)  to  (1),  we  see  that  the  vortieity  ( is  proportional  to  the  electron 
density  n by  the  relation  £ — 47rccZ?"!  n.  So,  in  the  electron  plasma  experiments, 
measxiring  electron  density  is  equivalent  to  measuring  vortieity. 


CHAOS  TO  CRYSTAL 

Figure  2(a)  shows  the  evolution  of  electron  density  (vortieity)  in  an  experiment 
performed  at  the  University  of  California,  San  Diego  [3].  Initially,  the  density  is 
distributed  in  a spiral,  resembling  a hot  tungsten  filament  that  is  used  to  create 
the  electrons.  Subsequently,  the  density  evolves  into  a system  of  N clumps  (dark 
spots)  immersed  in  a diffuse  background  plasma.  Each  density  clump  corresponds 
to  an  intense  vortex,  spinning  counter-clockwise,  and  the  diffuse  background  plasma 
corresponds  to  a low  level  of  background  vortieity. 

To  appreciate  the  influence  of  background  vortieity  on  the  vortex  mot  ion,  let  us 
first  consider  the  vortex  dynamics  without  background  vortieity.  In  the  absence 
of  background  vortieity,  the  vortex  dynamics  is  Hamiltonian.  Ignoring  minor  wall 
effects,  the  Hamiltonian  is  given  by  [4] 

H = ln  Ka:‘  ~ xs)2  + ('Ui  ~ Vj)2)]  - (3) 

i<J 
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FIGURE  2.  Spontaneous  formation  of  a vortex  crystal  in  an  initially  turbulent  electron  plasma 
[3].  (a)  Time  evolution  of  2-averaged  electron  density  (vorticity).  Darker  shades  indicate  higher 
density,  (b)  Number  of  vortices  versus  time,  (c)  Selection  of  vortex  crystals  obtained  from  initial 
conditions  similar  to  that  in  (a). 


where  Tj  is  the  circulation  (area  integral  of  vorticity)  of  the  ith  vortex,  ( Xi,yi ) is  the 
position  of  the  ith  vortex  in  a Cartesian  coordinate  system,  and  i,  j = 1,2,  ...N. 
The  canonically  conjugate  coordinates  are  ( qt,Pi ) = (v/T*  £*,  \/T7  y% ),  and  the 
dynamics  is  given  by  Hamilton’s  equations:  qt  = dPiH,  pl  = —dq,H. 

Because  the  vortex  dynamics  is  Hamiltonian  (without  background  vorticity),  it 
is  generally  chaotic.  During  this  chaotic  motion,  two  vortices  occasionally  become 
close  enough  to  merge.  Numerical  simulations  of  2D  Hamiltonian  vortex  dynamics, 
punctuated  by  occasional  mergers,  show  that  the  number  of  vortices  Nv  decays  as  a 
power  law  in  time,  t~*,  where  £ is  a positive  constant  of  order  unity  [5].  Figure  2(b) 
compares  the  expected  power-law  decay  to  an  electron  plasma  experiment.  The 
scatter  in  the  data  is  due  to  the  destructive  imaging  technique.  Initially,  Nv  exhibits 
the  expected  power-law  decay;  however,  when  t ~ 10ms,  mergers  cease. 

The  arrest  of  vortex  mergers  is  due  to  the  formation  of  a “vortex  crystal”  equi- 
librium [Fig  2(a),  100ms;  Fig.  2(c)].  A vortex  crystal  is  an  array  of  intense  vortices 
that  rotates  rigidly  in  a diffuse  background.  Jin  and  Dubin  proposed  that  vortex 
crystals  form  due  to  the  ergodic  mixing  of  background  vorticity  by  the  intense  vor- 
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tices  [6].  Using  a version  of  Lynden-Bell  statistics  [7],  they  showed  that  a vortex 
crystal  equilibrium  is  the  most  probable  state  that  can  result  from  such  mixing,  for  a 
given  energy,  angular  momentum  and  total  circulation  of  the  flow.  These  maximum 
entropy  states  are  compared  to  the  experiments  in  Ref.  [6] : they  accurately  predict 
the  arrangement  of  vortices,  and  details  of  the  background  vorticity  distribution.2 
Since  the  electron  plasma  experiments,  the  formation  of  vortex  crystals  has  been 
observed  in  various  numerical  simulations  of  2D  Euler  flow  [8,9]. 


MOTION  TO  THE  EXTREMA 

The  spontaneous  formation  of  vortex  crystals  in  2D  turbulence  demonst  rates  the 
importance  of  background  vorticity  in  spectacular  fashion.  Of  course,  other  effects 
of  background  vorticity  have  been  recognized  for  decades.  For  example,  it  is  well- 
known  that  vortices  move  toward  extrema  of  the  background  vorticity  distribution. 
The  details  of  this  process  have  been  studied  rigorously  in  dynamical  meteorology 
[12],  for  the  main  purpose  of  understanding  hurricane  trajectories.3  However,  a 
new  theory  was  required  to  explain  the  vortex  motion  in  electron  plasmas  [13]. 

Unlike,  terrestrial  hurricanes,  vortices  in  electron  plasmas  are  typically  exposed 
to  intense  background  shear.  In  this  sense,  electron  vortices  arc  more  like  Jovian 
hurricanes,  which  are  contained  in  strongly  sheared  zonal  jets  [10],  or  mesovor- 
tices  within  hurricanes  [9,11].  Intense  background  shear  creates  a disparity  in  the 
rates  at  which  vortices  of  opposite  sign  migrate  to  extrema  of  background  vorticity. 
Furthermore,  sufficiently  large  shear  will  arrest  the  migration.  These  results  are 
discussed  below. 

Figure  3(a)  is  a numerical  simulation  that  illustrates  vortex  motion  through  ax- 
isymmetric  background  vorticity,  similar  to  what  might  exist  in  an  electron  plasma. 
At  t = 0,  a clump  (black  dot)  and  a hole  (white  dot)  are  placed  at  the  same  radius 
tv  The  clump  is  a positive  vortex,  spinning  counter-clockwise,  whereas  the  hole  is 
a negative  vortex,  spinning  clockwise.  In  time,  both  vortices  redistribute  the  local 
background  vorticity.  In  response,  the  clump  moves  radially  inward,  to  the  peak  of 
background  vorticity,  whereas  the  hole  moves  radially  outward,  toward  a minimum 
of  background  vorticity. 

The  opposite  motion  of  clumps  and  holes  is  easily  understood.  Suppose  that 
there  is  only  one  vortex  of  circulation  F„,  and  that  the  vortex  is  point-like.  The 
vorticity  distribution  is  a sum  of  the  background  contribution  (b)  and  the  vortex 
contribution  (w);  i.e.,  C = (b  + r„<S(r  — f„).  It  is  well-known  that  the  2D  Euler 
equations,  with  rotationally  symmetric  boundary  conditions,  conserve  the  canonical 
angular  momentum, 


2*  The  final  number  of  vortices  and  their  individual  vorticity  profiles  are  not  predicted.  They  are 
taken  from  the  observed  vortex  crystals  and  fixed  in  computing  the  maximum  entropy  states. 

3)  The  background  vorticity  affecting  hurricane  motion  includes  planetary  vorticity,  as  well  as 
the  vorticity  associated  with  the  environmental  wind. 
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(a) 


FIGURE  3.  Vortex  motion  in  an  axisyrnmetric  background  vorticity  distribution,  (a)  Numerical 
simulation.  T = t x |S[r,j(0)]|.  (b,c)  The  redistribution  of  local  background  vorticity  by  a vortex 
increases  < r2  >b.  By  conservation  of  P$,  clumps  and  holes  must  react  oppositely. 

Pg  = J d2rC,r2  = Tb  < r2  >b  + r2v.  (4) 

Here,  Tp  > 0 is  the  total  circulation  of  the  background  flow,  and  < r2  >b  is 
the  (^-weighted  spatial  average  of  r2.  As  illustrated  in  Figs.  3(b)  and  3(c),  the 
redistribution  of  background  vorticity  by  the  vortex  tends  to  create  a local  plateau 
in  the  0-average  of  (V  This  increases  < r2  >b.  To  conserve  Pg,  a clump  (r„  > 0) 
must  decrease  rv  and  move  up  the  background  vorticity  gradient,  whereas  a hole 
(T„  < 0)  must  increase  rv  and  move  down  the  gradient. 

Let  us  now  consider  the  vortex-background  interaction  in  greater  detail.  The  Eu- 
ler equation  for  the  evolution  of  vorticity  can  be  written  as  two  separate  equations. 
The  first,  equation  is  for  the  vortex  motion, 


— f„  = vb(rv,0v,t). 


(5) 


Here,  vb  is  the  background  velocity  field;  that  is,  vb  = — x z,  where  V2^  = (j,. 
Since  the  background  vorticity  is  initially  axisyrnmetric,  i.e.,  C& (r,  t = 0)  = 0,(r), 
the  background  velocity  is  initially  circular,  i.e.,  vb(r,t  = 0)  = vb(r)0.  The  second 
equation  is  for  the  advection  of  background  vorticity, 

-QtCb  + (vb  + Vv)  ■ V(6  = o,  (6) 

where  vv  = z x r„(27r)-1V  ln(|r  — rv\)  plus  a small  correction  due  to  the  wall.  The 
vortex  trajectory  is  obtained  by  integrating  Eqs.  (5)  and  (6)  simultaneously. 
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FIGURE  4.  (Left)  Stream  lines  about  a retrograde  clump  and  a prograde  hole.  (Right)  The 
radial  speeds  of  retrograde  clumps  and  prograde  holes  of  various  strengths  (//;•„)  in  axisymrnetric 
background  vorticity.  Retrograde  clumps:  the  xs  and  diamonds  are  numerical  simulations,  the 
triangles  are  electron  plasma  experiments,  and  Eq.  (8)  is  computed  with  c = 0.43.  Prograde  holes: 
the  circles  are  numerical  simulations,  and  the  squares  are  electron  plasma  experiments. 


Neglecting  changes  to  Q [Eq.  (6)],  a single  vortex  merely  orbits  the  center  of 
background  vorticity,  with  angular  frequency  dOv/dt  = r~lvi,(rv).  However,  as  ex- 
plained previously,  the  background  vorticity  perturbation  causes  the  vortex  to  drift 
radially.  The  radial  velocity,  drv/dt , is  sensitive  to  whether  the  vortex  is  retro- 
grade or  prograde.  A vortex  is  ret  rograde  if  it  rotates  against  the  local  background 
shear,  whereas  a vortex  is  prograde  if  it  rotates  with  the  local  background  shear. 
Precisely,  let  S(r)  denote  the  local  shear-rate,  r d(vt,r~l)/dr.  Then,  r„/S(r„)  < 0 
for  a retrograde  vortex,  whereas  rv/S(rv)  > 0 for  a prograde  vortex.  Here,  we 
consider  examples  in  which  the  background  vorticity  decreases  monotonically  with 
r,  so  that  S(r)  is  negative.  Therefore,  in  these  examples,  clumps  arc  retrograde 
and  holes  are  prograde. 

The  stream  lines  about  a retrograde  clump  and  a prograde  hole  are  shown  in 
Fig.  4 (left).  There  are  two  stagnation  points  at  radial  distances  ±1  from  the  center 
of  the  retrograde  vortex,  where 


In  contrast,  there  are  no  stagnation  points  in  the  vicinity  of  a prograde  vortex. 
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Schecter  and  Dubin  [13]  derived  analytical  expressions  for  the  radial  velocities  of 
retrograde  and  prograde  vortices,  for  the  regime  in  which  l/rv  « 1.  This  regime 
corresponds  to  a weak  vortex  in  strong  background  shear.  For  a retrograde  vortex, 

Jfv  = ± | /2  ln  (C  rv/1)  - (8) 

where  +/—  is  for  clumps/holes,  and  c is  a constant  of  order  unity  that  depends  on 
the  particular  form  of  0>(r)-  For  a prograde  vortex, 


— r = ± — ^(r  ) l 2 
dt  v 4tt  dr[  v) 


(9) 


In  both  cases,  the  radial  speed  increases  with  the  local  background  vorticity  gra- 
dient, and  decreases  as  the  local  shear-rate  intensifies.  Note  that  the  prefactor 
(1/ 47r)  in  Eq.  (9)  is  much  less  than  the  prefactor  (7t/2)  in  Eq.  (8).  Moreover,  the 
logarithmic  factor  does  not  appear  in  Eq.  (9).  Thus,  as  l/rv  approaches  zero,  a 
prograde  vortex  will  move  orders  of  magnitude  slower  than  a retrograde  vortex  of 
equal  strength. 

The  plot  in  Fig.  4 shows  the  radial  vortex  velocity  ( drv/dt ) as  a function  of 
vortex  strength  ( l/rv ) for  a set  of  numerical  simulations,  and  for  several  electron 
plasma  experiments.  In  each  case,  a single  vortex  moved  through  a background 
vorticity  resembling  that  in  Fig.  3(a).4  The  simulations  are  described  in  Ref.  [13], 
and  the  experiments  are  described  in  Refs.  [14,15],  The  data  compares  favorably 
to  the  analytical  predictions  of  Schecter  and  Dubin.  Recently,  Kiwamoto  and  col- 
laborators used  a magnetized  electron  plasma  to  examine  the  motion  of  retrograde 
clumps  more  extensively  [16].  Their  data  also  agrees  with  Eq._(8). 

The  data  in  Fig.  4 are  for  cases  in  which  e(r)  = ~ (— ) /~  « 1 at  the  radial 

dr  \r  > dr 


position  rv  of  the  vortex.  However,  if  e(rv)  surpasses  a critical  value,  Eqs.  (8)  and 
(9)  become  invalid,  and  the  radial  velocity  of  the  vortex  drops  abruptly  to  zero.  For 
prograde  vortices,  Schecter  and  Dubin  estimated  that  the  critical  value  of  e(r„)  is 
unity.  For  retrograde  vortices,  Schecter  and  Dubin  estimated  that  the  critical  value 
of  e(rv)  is  of  order  \Z8ir2  ln(rv / /) . Both  estimates  are  consistent  with  numerical 
simulations  [13]. 


DANCING  WITH  THE  DIOCOTRON  MODES 

The  theoretical  work  of  Schecter  and  Dubin  [13]  neglects  the  interaction  of  the 
vortex  with  the  global  modes  of  the  background  vorticity,  i.e.,  the  diocotron  modes. 
While  usually  valid  (Fig.  4),  such  neglect  is  unjustified  if  the  orbital  frequency  of  a 
vortex  resonates  with  a neutral,  or  weakly  damped  diocotron  mode.  In  such  a case, 

b There  is  one  exception:  the  datum  with  error  bars  was  obtained  from  an  electron  plasma 
experiment  in  which  there  were  two  prograde  holes  [15]. 
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(a) 


FIGURE  5.  Example  in  which  the  radial  motion  of  a dump  is  controlled  by  diocotron  modes, 
as  opposed  to  a local  background  vorticity  perturbation,  (a)  Initial  vorticity  distribution,  (b) 
Diocotron  mode  amplitudes  and  radial  position  of  the  clump  versus  time,  in  units  of  27r/?-/IV  A 
sequence  of  resonances  causes  the  clump  to  merge  with  the  core  of  the  background  vorticity  [17]. 


there  can  be  significant  radial  motion  even  if  -r^(r„)  is  zero,  or  if  s(rv)  is  above  its 

dr 

critical  value.  The  following  example  illustrates  this  point. 

Consider  the  initial  vorticity  distribution  in  Fig.  5(a).  The  background  vorticity 
0,(r)  consists  of  a uniform  core  (r  < R),  and  an  outer  skirt  (r  > /?),  where  the 
vortex  resides.  By  assumption,  the  vorticity  in  the  outer  skirt  is  negligible,  so  that 
r b ~ 7ri?2Cfc(0).  In  addition,  dQ/dr  ~ -C/>(0)A(r  - R). 

The  background  vorticity  perturbation  will  be  concentrated  at  the  core  radius 
R.  Expanding  this  perturbation  in  a Fourier  series  yields 


C b - C b(r)  + - R)  cos [m6-<pm(t.)], 


(10) 


where  Yb/itR  is  a convenient  normalization  factor.  The  Fourier  components  are 
the  “diocotron  modes”  of  the  background  vorticity.  Exciting  a diocotron  mode 
corresponds  to  displacing  (rn  — 1)  or  deforming  (m  > 2)  the  core.  In  the  following, 
it  will  be  assumed  that  the  vortex  is  weak  and  that  the  mode  amplitudes  are  small; 
i.e.,  7 = Tu/rj  <<  1,  and  n,„  <<  1.  Changes  to  the  axisymmetric  component  of 
background  vorticity  are  of  order  a2,,  and  are  neglected  in  Eq.  (10). 

The  vortex-mode  dynamics  may  be  put  in  Hamiltonian  form  [17].  The  canoni- 
cally conjugate  variables  for  the  mlh  mode  are  (r/m,  p,„ ) = (—(/?,„,  u'lJrn),  and  the 
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canonically  conjugate  variables  for  the  vortex  are  (qv,  pv)  = {—yfy  0V,  y/j  rl/R2). 
The  equations  of  motion  are  Hamilton’s  equations:  (qm,  pm)  = ( dVmH , — dqmH), 
and  (qv,  pv)  = (dPvH,  —dqvH),  where  (ignoring  wall  effects) 


H 2tt/?2  X 

I OO 

< 7 ln(p„)  + ]P(m  - 1)  p„ 

l 7715=1 


(11) 


The  first  term  in  braces  accounts  for  the  advection  of  the  vortex  in  the  unperturbed 
circular  flow.  The  middle  term  accounts  for  the  free  propagation  of  diocotron 
modes,  in  the  absence  of  the  vortex.  The  final  term  couples  the  vortex  motion  to 
the  modes.  The  Hamiltonian  in  Eq.  (11)  neglects  mode-mode  interaction  terms, 
which  are  of  order  afn. 

Of  course,  the  vortex-mode  interaction  conserves  canonical  angular  momentum. 
Minus  the  constant  unperturbed  background  contribution,  the  (dimensionless) 
canonical  angular  momentum  can  be  written  as  follows: 


oo  2 

Po  = + 7^-  (12) 

771=1 

Suppose  that  a clump  (7  > 0)  excites  a mode.  Then,  increases.  To  conserve 
Pg,  rv  must  decrease.  In  contrast,  if  a hole  (7  < 0)  excites  a mode,  rv  must  increase 
to  conserve  Pg. 

Figure  5(b)  illustrates  the  vortex-mode  dynamics  governed  by  the  Hamiltonian 
in  Eq.  (11).  At  t = 0,  a clump  (7  = .0025)  is  placed  at  rv  — r0  = 1.418/?,  and 
the  mode  amplitudes  are  all  zero.  The  clump  principally  orbits  the  center  of  the 
background  with  angular  velocity  ddv/dt.  ~ vb(rv)/rv  = Fb/2-nrl.  At  rv  = 1.418.R, 
this  angular  velocity  equals  the  natural  phase  velocity  of  the  m — 2 diocotron  mode. 
Consequently,  the  clump  resonantly  excites  the  m — 2 mode.  By  conservation 
of  Pg,  this  resonance  causes  the  clump  to  move  radially  inward.  As  the  clump 
approaches  the  core,  d0v/dt  increases.  This  causes  resonances  with  higher  order 
diocotron  modes:  first  the  m = 3,  then  the  m = 4,  etc.  The  sequence  of  resonances 
incrementally  pulls  the  clump  into  the  core. 

Lansky,  O’Neil  and  Schecter  [17]  derived  a condition  for  which  the  clump  will 
merge  with  the  core.  For  the  case  in  which  the  initial  clump  position  is  r„  = 
1.418/?,  7 must  be  greater  than  about  0.002.  Otherwise,  the  initial  oscillations  of 
rv  [Fig.  5(b)]  will  not  bring  the  clump  close  enough  to  the  core  to  resonantly  excite 
higher  order  modes.  Electron  plasma  experiments  have  yet  to  verify  this  merger 
condition. 

However,  Durkin  and  Fajans  [18]  have  used  a magnetized  electron  plasma  to 
examine  the  interaction  of  diocotron  modes  with  a clump  that  is  already  inside  the 
core  (r„  < /?).  In  this  case,  there  are  no  vortex-mode  resonances.  Nevertheless, 
the  clump  can  weakly  excite  the  diocotron  modes.  This  creates  a wave  along  the 
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edge  of  the  core.  Eventually,  the  wave  breaks.  The  wave-breaking  process  alters 
the  background  vorticity  profile,  and  can  generate  one  or  more  vorticity  holes. 
The  measured  wave-breaking  time  is  proportional  to  7_1ln(7_1),  provided  that 
r„  < 0.7 R.  This  result  agrees  with  a recent  theory  by  Jin  and  Dubin  [19]. 

Clearly,  background  vorticity  enriches  2D  vortex  dynamics.  Background  vorticity 
can  cause  the  spontaneous  formation  of  vortex  crystals  (Fig.  2),  and  can  spat  ially 
separate  vortices  of  opposite  sign  (Fig.  3).  Recent  work  in  non-neutral  plasma 
physics  has  shed  some  light  on  2D  vortex  dynamics  with  background  vorticity,  but 
surely  more  surprises  await. 
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